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I.  INTRODUCTION 


Terminal  ballist icians  have  the  requirement  for  knowing  the  mass  and 
velocity  of  fragments  generated  by  munition  impacts  on  armor  targets.  In  the 
past,  a  great  deal  of  human  effort  has  been  expended  in  order  to  acquire  this 
information  and  often  the  accuracy  of  the  the  data  is  poor.  For  example,  the 
fragments  produced  behind  target  plates  have  been  recovered  by  stopping  them 
in  celotex  and  the  masses  of  the  fragments  and  their  direction  of  travel  can 
be  determined  quite  accurately  after  a  lot  of  tedious  work,  but  the  speed  of 
the  fragments  can  only  be  estimated. 

An  alternative  method  to  capturing  the  fragments  in  celotex  is  to 
capture  images  of  the  fragments  in  flight  on  film  by  pulsing  orthogonal  X-ray 
tubes  at  known  times.  Again,  a  great  deal  of  human  effort  is  required  to 
analyze  the  film,  often  with  unsatisfactory  results.  The  problem  is  that  it  is 
difficult  to  correlate  orthogonal  images  of  fragments  (required  in  order  to 
determine  fragment  mass)  and  even  more  difficult  to  identify  multiple  fragment 
images  in  the  same  view  produced  by  a  multiflash  sequence  which  is  needed  to 
obtain  accurate  fragment  velocities. 

The  difficulties  associated  with  fragment  data  acquisition  and  analysis 
can  now  be  alleviated  by  recent  advancements  in  computer  and  video  technology. 
It  is  now  feasible  to  automate  the  data  acquisition  and  analysis  system  so 
that  the  amount  of  human  effort  will  be  greatly  reduced  and  the  accuracy  and 
speed  in  the  analysis  substantially  enhanced.  This  report  describes  recent 
advancements  in  fragment  image  analysis  and  presents  the  technical  basis  for  a 
proposed  Automatic  Data  Acquisition  System  (ADAS). 


II.  ANALYSIS  OF  FRAGMENT  IMAGES 

The  mathematical  procedure  for  analyzing  the  fragment  image  data  in  the 
ADAS  will  be  essentially  the  same  as  that  used  for  the  film  image  data. 
However,  in  the  ADAS,  the  images  will  be  recorded  electronically  and  the 
analysis  performed  automatically.  In  order  to  show  the  feasibility  of  the 
ADAS,  it  is  first  necessary  to  discuss  the  procedure  for  calculating  the  mass 
and  velocity  of  the  fragments  from  the  film  data.  The  basis  for  this  analyti¬ 
cal  procedure  has  been  reported  previously, ^  but  some  modifications 
which  constitute  improvements  in  the  technique  are  presented  in  this  section. 

*A.  L.  Arbuckle ,  E.  L.  Herr,  A.  J.  Ricchiazzi,  "A  Computerized  Method  of 
Obtaining  Behind- the-Tar get  Data  From  Orthogonal  Flash  Radiographs,  rr  BRL 
MR  2264,^ January  1973  (AD  908362L). 
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II. 1  SCALING  THE  FILM  COORDINATE  POINTS  TO  THE  SPATIAL  COORDINATES 


p 

The  following  analysis  is  based  on  the  typical  X-ray  system  used  at  the 
BRL.  The  geometric  relationship  of  the  orthogonal  X-ray  tubes  to  the  film  is 
depicted  in  Figure  1.  As  shown,  the  origin  of  the  coordinate  system  is  defined 
to  be  the  intersection  of  the  two  film  planes  (the  X-Z  plane  and  Y-Z  plane) 
with  the  X-Y  plane.  The  X-Y  plane  is  formed  by  the  centerlines  of  the  orthog¬ 
onal  X-ray  tubes.  The  centerline  of  the  horizontal  tube  is  oriented  parallel 
to  the  X  axis  with  the  tube  facing  the  Y-Z  plane  (vertical  film  plane)  and  is 
a  distance  Y^  above  the  X  axis.  The  forward  tip  of  the  horizontal  tube  is  a 
distance  X^  from  the  Y-Z  plane.  The  centerline  of  the  vertical  tube  is 
oriented  parallel  to  the  Y  axis  with  the  tube  facing  the  X-Z  plane  (horizontal 
film  plane)  and  is  a  distance  Xff  away  from  the  Y  axis.  The  forward  tip  of  the 
vertical  tube  is  a  distance  Yhf  from  the  X-Z  film  plane. 


Reference  marks  needed  to  analyze  the  image  data  are  made  on  the  film  by 
fiducial  wires.  A  fiducial  wire  parallel  to  the  Z  axis  is  located  directly  in 
front  of  the  vertical  film  at  a  distance  Yff  above  the  Z  axis.  Another  fidu¬ 
cial  wire  parallel  to  the  Z  axis  lies  directly  above  the  horizontal  film  at  a 
distance  X^  away  from  the  Z  axis.  Fiducial  wires  (not  visible  in  Figure  1) 
are  also  strung  along  the  X  axis  and  the  Y  axis.  The  point  of  intersection  of 
each  pair  of  fiducial  wires  lies  on  the  centerline  axis  of  the  orthogonal  X- 
ray  tube  which  faces  that  point. 


In  order  to  determine  the  spatial  coordinates  (coordinates  of  points 
lying  in  the  space  directly  in  front  of  both  film  planes)  of  certain  points  on 
the  fragments  from  corresponding  points  on  their  respective  images  on  the 
orthogonal  films,  it  is  necessary  to  scale  the  film  coordinates  of  the  points 
on  the  film.  (The  coordinates  of  points  located  on  the  film  surface  will  be 
referred  to  as  film  coordinates  and  have  the  same  reference  origin  as  the 
spatial  coordinates.)  The  reason  that  the  film  coordinates  must  be  scaled  to 
the  spatial  coordinates  is  due  to  the  fact  that  each  X-ray  tube  acts  as  a 
point  source  of  radiation.  Therefore,  any  point  in  space  between  the  tube  head 
and  the  film  except  those  points  lying  on  the  X-ray  tube  centerline  will  be 
nonisometrically  projected  onto  the  film. (The  points  that  lie  on  the 
centerline  are  isometrically  projected.)  The  derivation  of  the  equations  for 
evaluating  the  appropriate  scaling  factors  to  translate  the  film  coordinates 

^ Chester  Grabarek ,  Louis  Herr,  "X-Ray  Multiflash  System  for  Measurement  of 
Projectile  Performance  at  the  Target,"  BRL  TN  1634,  September  1966  (AD  807619). 
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Geometric  Relationship  Of  Orthogonal  X-Ray  Tubes  To  The  Film  Planes 


to  the  spatial  coordinates  is  described  in  the  remainder  of  this  section. 


The  equations  from  which  the  scaling  factors  are  derived  result  from  the 
geometric  relationship  of  the  X-ray  tubes  to  the  film  planes.  Figure  2  is  an 
expanded  view  of  Figure  1  in  the  region  close  to  the  origin.  Not  shown  is  the 
point  X  ,Y_,Zp  in  space  which  lies  at  the  intersection  of  the  lines  drawn  from 
each  X-rayP tube  to  its  corresponding  film  image  point.  The  coordinates  of  the 
vertical  film  image  point  are  YVZW  and  the  coordinates  of  the  horizontal  film 
image  point  are  X1fZh.  The  following  equations  are  derived  based  on  the  rules 
for  similar  triangles: 


Y1 


(1) 


and 


(2) 


Two  more  equations  can  be  derived  in  a  similar  manner  based  on  the  geometry 
shown  in  Figure  3.  These  equations  are: 


xhf 


Yi  -  y 


ff 


(3) 


and  Yhf 

xff  -  X1 


(4) 


From  these  equations,  a  scaling  factor  for  each  film  plane  can  be 
derived.  The  scaling  factor  is  defined  as  the  value  needed  to  scale  the  film  Z 
coordinate  value  to  the  spatial  Z  coordinate  value.  Adopting  the  symbol  used 
in  Reference  1  for  the  scaling  factor,  the  scaling  factor  for  the  vertical 
film  plane  is  by  definition: 


K 


v 


zp/zv 


(5) 


The  equation  for  calculating  the  value  for  Ky  is  found  by  solving  Equations  1, 
3  and  4  simultaneously.  This  is  done  by  taking  Equation  3,  solving  for  Xp  and 
substituting  the  result  in  Equation  4.  The  expression  found  by  cross 
multiplying,  expanding  and  solving  for  Yp  is  then  substituted  into  Equation  1 
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Figure  2  Expanded  View  Of  The  Geometry  Close  To  The  Origin 


VERTICRL 
X-RRY  TUBE 


Ul 

a:  m 

ZD 
~Z-  h- 


12 


L _ 


Figure  3  View  Of  The  X-Y  Plane  With  Z=0 


and  rearranged  in  accordance  with  Equation  5.  The  result  is: 


Yhf^Xhf  ~  X-P  ~  Yff^Xff  ~  X1^  ^ 

XhfYhf  +  (Xff  *  X1)(Y1  *  Yff^  ' 


Similarly  to  Equation  5,  the  scaling  factor  for  the  horizontal  film  is  bv 
definition : 


Kh  =  Zp/Zh  .  (7) 

Following  a  similar  nrocedu^e  as  that  used  to  obtain  Equation  b,  hut  using 
Equations  ?,  7,  4  and  7,  will  yield  the  following  equation. 


x  _  xhf(Yhf  -  Yi)  +  Xff-Yi  ~  Yff) 
^  "  XhfYhf  +  (Xff  "  X1)(Y1  "  Yff} 


(8) 


The  equations  for  scaling  the  film  coordinate  points  to  the  sDatial 
coordinate  Doints  using,  both  K  factors  are  derived  as  follows.  Rearranging 
Equations  E  and  7  results  in  the  equations  for  calculating  7,p  which  are: 


v  v 


(9a) 


or 


=  Xh-h 


(Qb) 


That  is,  Zp  can  be  calculated  from  the  Z  film  coordinate  value  from  either 
orthogonal  film  when  using  the  appropriate  K  factor.  By  substituting  Equation 
9a  into  Equation  1  and  rearranging,  the  equation  to  calculate  the  value  for  Yp 
is : 


YP  =  Yff  +  MYi  -  Yff>  >  (10) 

and,  similarly,  by  substituting  Equation  9b  into  Equation  ?  and  rea^rangi ng, 
the  equation  for  calculating  the  value  of  Xp  is: 

Xp  =  Xff  -  Kh(Xff  -  X-,)  .  (11) 


In  general,  the  values  for  Ky  and  Kh  are  not  equal  but  there  are  particu¬ 
lar  spatial  points  for  which  they  are  equal.  These  points  can  be  determined  by 
setting  Equation  6  equal  to  Equation  8  and  solving  for  X-|  for  a  given  Y^. 
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To  obtain  maximum  accuracy  in  evaluating  the  X  and  the  Y  spatial  coordi¬ 
nate,  the  respective  K  factor  should  be  used.  However,  there  are  some  situa¬ 
tions  in  which  this  can  not  be  done.  For  example,  when  the  target  is  at  an 
oblique  angle  to  the  vertical  plane,  any  fragments  in  the  region  directly  in 
front  and  the  region  directly  behind  the  target  will  be  masked  by  the  image  of 
the  target  plate  on  the  horizontal  film.  For  these  fragments,  the  only  images 
appearing  will  be  on  the  vertical  film.  In  such  cases,  an  assumption  must  be 
made  about  determining  the  spatial  locations  of  these  fragments  since  there  is 
insufficient  information  to  calculate  their  locations  directly.  The  assumption 
usually  made  is  that  the  fragments  are  in  the  same  vertical  plane  as  the  shot 
line  (the  penetrator  trajectory  before  impacting  the  target).  Therefore,  the 
scaling  factor  (Ky)  that  was  evaluated  to  obtain  the  striking  velocity  is  used 
to  scale  the  coordinates. 


II. 2  MATCHING  ORTHOGONAL  FILM  IMAGES  AND  CALCULATING  FRAGMENT  MASS 


Before  the  spatial  location  of  a  specific  fragment  can  be  determined  from 
the  images  in  the  orthogonal  films,  it  is  necessary  to  determine  which  image 
in  each  film  was  caused  by  that  fragment.  A  pair  of  orthogonal  images  that  can 
be  assigned  to  an  individual  fragment  is  called  a  match.  The  procedure  for 
matching  images  is  based  on  the  fact  that  both  film  planes  share  a  common  axis 
-  the  Z  axis.  Therefore,  the  extreme  endpoints  in  the  Z  direction  of  the 
fragment  in  space  appear  as  the  extreme  endpoints  in  the  Z  direction  of  the 
image  on  each  orthogonal  film.  Maximum  accuracy  in  evaluating  the  K  factors 
(and,  therefore, the  spatial  coordinates)  occurs  for  those  images  where  there 
is  a  unique  Y  value  and  a  unique  X  value  on  the  respective  films  to  pair  to 
each  Z  endpoint.  When,  for  example,  there  are  a  range  of  Y  values  for  an  image 
in  the  vertical  film  that  can  be  paired  to  a  Z  endpoint,  the  values  calculated 
for  the  K  factors  will  vary  depending  on  what  value  is  selected  over  the  range 
of  Y  values  but  the  values  for  the  K  factors  will  vary  only  slightly.  The 
reason  for  this  is  that  the  magnitude  of  the  change  in  the  Y  values  is  small 
compared  to  the  head-to-film  distance.  For  example,  assuming  that  Yhf=Xhf=51J, 
Yff=Xff=8,  X i  =  7  and  Yi  =7,  then  Kv=0.8734l.  Allowing  Y i  to  vary  over  a  range  of 
1  (which  represents  an  unusually  large  range)  so  that  Y i  =  6 ,  then  Ky=0. 89224 
or  a  change  of  about  2%.  This  would  mean  a  possible  2 %  error  in  determining 
the  spatial  Z  coordinate  value. 

The  first  step  in  the  process  for  matching  orthogonal  images  is  the 
selection  of  the  minimum  Z  endpoint  of  each  pair  of  orthogonal  images,  comput¬ 
ing  the  K  factor  for  each  endpoint  and  then  scaling  the  7,  values  to  the 
spatial  coordinates  by  using  Equations  9a  and  9b.  If  these  scaled  Z  values  are 
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within  some  tolerance  value  (typically  10  mm)  -  to  make  allowance  for  slight 
errors  in  measuring  the  coordinate  values  -  then  the  image  endpoints  were 
determined  by  the  single  spatial  point  located  at  the  minimum  Z  coordinate 
value  of  the  fragment.  The  next  step  is  to  process  the  maximum  Z  endpoints  of 
the  images  in  an  identical  manner.  When  both  the  minimum  Z  endpoints  and  the 
maximum  Z  endpoints  are  within  the  tolerance  value,  the  matching  process  for 
this  pair  of  images  can  be  terminated  if  all  of  the  film  images  are  well 
separated.  Usually,  the  images  are  not  well  separated,  so  two  more  steps  are 
taken  to  make  a  final  selection  on  which  images  will  be  declared  a  match. 

The  next  step  is  to  test  to  see  if  the  sign  of  the  difference  between  the 
two  minimum  Z  endpoints  is  the  same  sign  as  the  difference  between  the  two 
maximum  Z  endpoints  (when  taken  in  the  same  order).  This  is  based  on  the 
theory  that  the  error  in  measuring  the  Z  endpoints  of  each  image  will  have  the 
same  sign  -  that  is,  be  displaced  in  the  same  direction.  This  error  can  occur 
when  the  fiducial  wires  are  shifted  out  of  position  when  inserting  the  film 
cassette  and  are  not  returned  to  their  normal  position.  This  is  manifested  by 
the  fiducial  wire  images  not  being  perpendicular  to  each  other. 


The  final  step  in  the  matching  process  is  to  select  the  pair  of  images 
for  which  the  Z  length  of  each  image  is  most  nearly  alike.  The  Z  length  is  the 
maximum  Z  endpoint  minus  the  minimum  Z  endpoint  where  the  Z  endpoints  are  the 
values  obtained  from  Equations  9a  and  9b. 


To  summarize,  the  matching  procedure  is  outlined  as  follows: 

(1)  Select  an  image  from  the  vertical  film  (List  1). 

(2)  Select  an  image  from  the  horizontal  film  (List  2). 

(3)  Compute  the  K  factors  for  the  minimum  Z  endpoints  of  the  two  images. 

(4)  Adjust  the  Z  values  using  the  appropriate  K  factor. 

(5)  Take  the  absolute  value  of  the  difference  between  the  two  adjusted 

Z  values. 

(6)  If  this  difference  is  greater  than  the  selected  tolerance  value 

(typically  10  mm.)  then  return  to  Step  2  selecting  a  new  image 
from  List  2  or  if  it  is  depleted  then  Step  1 . 


15 


(7)  Repeat  Steps  3  through  6  using  the  maximum  Z  endpoints. 


(8)  For  those  that  meet  the  tolerance  test  for  both  minimum  and  maximum 

endpoints,  then  test  to  see  that  the  displacement  between  the  two 
minimum  Z  points  is  in  the  same  direction  as  the  displacement  in  the 
two  maximum  endpoints.  If  not,  then  return  to  Step  2  or,  if  appro¬ 
priate,  Step  1. 

(9)  More  than  one  image  may  qualify  as  a  match  based  on  the  above  tests. 

The  final  match  is  determined  by  selecting  the  pair  of  images  that 
are  most  nearly  identical  in  Z  length. 

For  each  pair  of  matched  images,  an  estimate  of  the  volume  of  the 
fragment  can  be  determined.  Briefly,  the  volume  is  estimated  by  evaluating  the 
area  of  the  matched  image  in  each  orthogonal  film,  constructing  a  rectangle  of 
equivalent  area  as  that  of  the  image  and  then  scaling  both  rectangles  to  the 
spatial  position  of  the  fragment.  This  results  in  a  rectangular  parallelepiped 
representing  the  maximum  possible  volume  of  the  fragment.  Since  this  is  the 
maximum  possible  volume,  an  empirically  determined  percentage  of  the  volume  is 
taken  to  represent  the  actual  volume  of  the  fragment. 

The  first  step  in  determining  the  area  of  an  image  is  to  outline  the 
image  with  a  series  of  connected  straight  line  segments  as  depicted  in  Figure 
M.  Then, proceeding  in  a  clockwise  direction  around  the  boundary,  the  end 
points  of  the  straight  line  segments  are  digitized  (the  coordinates  of  the 
endpoints  are  evaluated)  and  the  values  substituted  in  one  of  the  following 
equations,  depending  on  which  film  plane  is  being  processed.  This  procedure  is 
continued  until  the  outline  of  the  image  has  been  traced  back  to  the  starting 
point. 

N 

Ayz  =  H  E  (Y1ZU1  -  li.lV  ('2) 

i  =  l 

N 

or  Axz  =  ^  £cXiZi+1  -  X±^Z±)  (13) 

i  =  l 

The  centroid  of  the  image  is  also  evaluated  since  it  is  used  to  position 
the  rectangle  which  will  represent  the  image.  The  centroid  for  an  image  in  the 
vertical  film  is  calculated  using  the  following  equations: 

-1  N 
Yo  =  i- 

y  i  =  l 


16 


and 


1 


(15) 


Jcv 


6A 


N 

E 


(Y, 


i+1 


Y,)(Z t  + 


ZiZi+1 


zLi) 


yz  i=l 


The  centroid  for  an  image  in  the  horizontal  film  is  calculated  using  the 
following  equations: 
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(17) 


The  signs  involved  in  the  equations  for  the  centroid  require  a  clockwise 
rotation  as  the  boundary  line  segments  are  being  traced.  The  summations  invol¬ 
ved  in  the  centroid  equations  are  evaluated  at  the  same  time  that  the  summa¬ 
tion  for  the  area  equation  is  computed.  Upon  returning  to  the  starting  point, 
the  value  for  the  area  will  have  been  determined  and  the  calculations  for  the 
centroid  can  be  completed.  For  Equations  12  -  17,  the  N+1  point  is  the  same  as 
the  starting  point. 


The  position  of  the  sides  of  the  rectangle  which  will  represent  the 
image  is  determined  by  aligning  the  sides  with  the  Z  endpoints.  The  "top"  and 
"bottom"  of  the  rectangle  is  determined  by  centering  the  rectangle  in  the  Y 
(or  X)  direction  on  the  centroid.  This  is  illustrated  in  Figure  5.  The  loca¬ 
tion  of  the  point  to  represent  the  "top"  of  the  rectangle  is  chosen  to  be  the 
point  aligned  "above"  the  Z  coordinate  of  the  centroid.  Likewise,  the  point 
aligned  "below"  the  Z  coordinate  of  the  centroid  represents  the  "bottom"  of 
the  rectangle.  The  equations  involved  in  establishing  these  points  for  the 
vertical  film  are  the  following: 


Az 

Ay 

"top"  point 

and  the  "bottom"  point 

Similar  equations  apply  for  the  image 


(7  -7  ) 

'max  min'  ’ 

(18) 

Ayz/AZ  , 

(19) 

(ycv  +  Ay/2,  zcv) 

(20) 

(Ycv  -  Ay/2,  zcv)  . 

(21) 

in  the  horizontal  film. 
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To  summarize,  there  are  five  points  associated  with  each  film  image. 
These  five  are  the  minimum  Z  end  point,  the  maximum  Z  endpoint,  the  centroid, 
the  "top  center"  point  and  the  "bottom  center"  point.  For  each  pair  of  matched 
images,  the  five  points  from  each  image  are  scaled  using  Equation  9a,  9b,  10 
or  11  appropriately  to  form  a  rectangular  parallelepiped  in  the  spatial  coor¬ 
dinate  system.  (The  K  factor  which  is  used  to  scale  the  "top"  point,  the 
centroid  and  the  "bottom"  point  is  the  average  value  of  the  K  factors  for  the 
Z  endpoints.)  The  parallelepiped  represents  the  maximum  possible  volume  of  the 
fragment.  To  represent  the  actual  volume  of  the  fragment,  a  fractional  part  of 
volume  of  the  rectangular  parallelepiped  is  required.  A  satisfactory  value  for 
this  percentage  value  will  be  determined  in  the  next  subsection. 

Once  the  volume  of  the  fragment  has  been  calculated,  the  mass  is  computed 
by  multiplying  the  volume  by  the  density  of  the  fragment,  (if  the  target  and 
penetrator  are  made  from  different  materials,  it  may  be  necessary  to  guess 
which  density  should  be  used.)  The  equation  to  compute  the  mass  of  the  frag¬ 
ment  is: 


Mf  =  q  p  AX'  AY*  A  Z»  .  (22) 

where:  q  is  the  fractional  volume  value,  is  the  fragment  density  and  the 
primes  indicate  that  the  dimensions  of  the  rectangular  parallelepiped  are 
determined  from  the  scaled  coordinate  points  (i.e.,  the  spatial  coordinate 
points) . 


II. 3  A  TEST  OF  THE  METHODOLOGY  FOR  MATCHING  IMAGES  AND  THE  DETERMINATION 

OF  THE  FRACTIONAL  VALUE  q 

To  demonstrate  how  well  the  rnethodolgy  for  matching  images  works,  six 
irregular  shaped  fragments  ranging  in  mass  from  9.75  grams  to  1M8.5  grams  were 
selected  from  an  assortment  of  fragments  recovered  from  a  firing  range.  These 
fragments  were  mounted  at  various  positions  on  a  piece  of  styrofoam  in  such  a 
way  that  the  projected  images  on  the  each  of  the  orthogonal  films  would  be 
completely  separated  from  each  other.  A  pair  of  orthogonal  X-ray  tubes  were 
pulsed  and  the  resulting  images  formed  on  the  two  orthogonal  films  were 
analyzed  using  the  procedure  described  in  the  previous  subsection.  The 
measured  Z  endpoints  and  the  calculated  coordinates  ("top,11  centroid  and 
"bottom")  representing  the  rectangle  with  an  equivalent  area  to  that  of  each 
image  on  the  film  are  tabulated  in  Tables  1  and  2  along  with  the  computed 
areas  of  the  images  on  the  film.  Table  3  shows  the  results  of  matching  the 
images.  Since  the  images  were  well  separated,  a  unique  match  was  found  for 
each  of  the  images. 
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TABLE  1  VERTICAL  FILM  COORDINATE  VALUES  FOR  A  TEST  CASE 
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TABLE  2  HORIZONTAL  FILM  COORDINATE  VALUES  FOR  A  TEST  CASE 
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TABLE  3  MATCHED  IMAGE  PAIRS  (SCALED  SYSTEM  COORDINATES  )  FOR  A  TEST  CASE 
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The  maximum  possible  mass  for  each  matched  image  pair,  found  by  multi¬ 
plying  the  density  of  steel  (7.78  g/cc)  by  the  volume  of  the  rectangular 
parallelepiped,  is  tabulated  in  the  fifth  column  of  Table  4.  The  ratio  of  the 
measured  mass  to  the  calculated  mass  is  tabulated  in  the  sixth  column  of  Table 
4.  The  average  value  of  these  ratios  is  0.47.  Therefore,  it  is  suggested  that 
the  mass  of  an  irregular  shaped  fragment  can  be  estimated  from  the  following 
equation. 


Mf  =  0.47  pf  AX'  AY'  AZ'  . 


(23) 


The  percentage  error  in  calculating  the  mass  using  this  equation  ranges 
from  3.6*  to  17.7*  for  these  six  fragments.  The  average  error  is  11.4*. 


II. 4  ESTIMATING  FRAGMENT  VELOCITY  BASED  ON  ORTHOGONAL  FILM  IMAGES 


Calculating  the  velocity  of  a  fragment  requires  establishing  the  position 
of  the  fragment  at  two  locations  and  the  time  interval  which  elapsed  while  the 
fragment  travelled  from  the  first  location  to  the  second.  Two  methods  for 
calculating  fragment  velocity  will  be  described. 


The  first  method  requires  only  one  X-ray  station  (one  set  of  orthogonal 
X-ray  tubes  )  and  makes  use  of  two  simplifying  assumptions.  The  first 
assumption  is  that  all  fragments  leave  from  a  single  point  and  the  second 
assumption  is  that  they  all  leave  from  that  point  simultaneously.  The 
following  discussion  concerning  the  magnitude  of  the  error  involved  by  making 
these  assumptions  will  be  based  on  the  situation  where  fragments  leave  the 
rear  surface  of  the  target  as  a  result  of  complete  penetration. 


The  magnitude  of  the  error  involved  in  making  the  first  assumption  (frag¬ 
ments  leaving  from  a  single  point)  is  small  since  the  dimensions  of  the  exit 
hole  are  small  compared  to  the  distance  the  fragments  have  travelled  by  the 
time  the  X-ray  tubes  are  flashed.  This  single  point  is  chosen  to  be  the  center 
of  the  exit  hole  in  the  original  plane  of  the  rear  surface  of  the  target.  It 
is  denoted  by  the  coordinates  Xg,  Yg,  Zg  where  the  reference  origin  is  the 
same  as  that  of  the  film  and  spatial  coordinate  system. 


TABLE  4  CALCULATED  MASS  FROM  FLASH  RADIOGRAPHS  FOR  A  TEST  CASE 
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The  magnitude  of  the  error  involved  in  making  the  second  assumption  is 
more  difficult  to  ascertain.  It  is  thought  that  the  period  of  time  required 
for  the  fragments  to  reach  their  maximum  velocity  is  small  compared  to  the 
interval  of  time  culminating  in  the  flash  of  the  X-ray  tubes. 

The  magnitude  of  the  total  error  is  dependent  on  the  method  used  to 
establish  the  beginning  of  the  time  interval  to  be  used  in  the  velocity 
calculation.  The  recommended  method  for  establishing  the  beginning  of  the  time 
interval  is  to  place  a  trigger  break  screen  approximately  one  inch  (2.54  cm) 
behind  the  target  parallel  to  the  rear  surface  by  using  a  layer  of  flexible 
foam  material.  The  time  interval  used  in  the  velocity  calculation  is  the  time 
between  the  instant  the  screen  is  broken  (which  triggers  a  time  counter)  to 
the  instant  the  X-ray  tubes  are  pulsed  (stopping  the  time  counter).  This  time 
interval  will  be  denoted  as  At  in  the  equation  below. 


Using  this  method,  the  equation  for  calculating  the  velocity  of  fragments 
which  have  been  identified  by  matched  image  pairs  is: 


vf  - 


v 


(Xp  -  Xe)2  *  (Ip  -  Ie)z  * 


«Zph  *  V)/2  - 


(24) 


At 


where  X^,  and  Yp,  ZpV  are  the  scaled  coordinates  (spatial  coordinates)  of 

the  centroids. 


The  second  method  for  computing  fragment  velocities  requires  two  orthogo¬ 
nal  X-ray  stations  (  a  pair  of  orthogonal  X-ray  tubes  at  each  station)  behind 
the  target.  Since  the  same  film  is  exposed  by  the  flash  from  each  X-ray 
station,  the  images  of  the  fragments  produced  by  the  flash  from  the  first  X- 
ray  station  lie  close  to  the  images  produced  by  the  flash  from  the  second  X- 
ray  station.  It  is  frequently  difficult  to  determine  to  which  station  each 
image  should  be  assigned.  (This  problem  will  not  arise  in  the  ADAS  because  the 
images  from  each  X-ray  station  will  be  separate  from  those  of  any  other 
station. ) 


The  orthogonal  images  are  matched  at  each  station  in  the  manner  described 
in  the  previous  subsection.  Then  the  images  at  one  station  must  be  matched  to 
the  images  at  the  other  station.  This  is  not  easy  to  do  since  the  image  of  a 
fragment  at  one  station  does  not  necessarily  look  the  same  as  its  related 
image  at  the  subsequent  station  due  to  rotation  and  tumbling.  One  way  to  make 
the  match  is  to  determine  the  trajectory  of  a  fragment  at  the  first  station 
based  on  a  virtual  reference  point  (one  point  for  all  the  fragments)  and  look 
for  a  fragment  at  the  second  station  which  lies  on  that  trajectory.  When  a 
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match  can  be  made,  the  velocity  can  be  determined  by  evaluating  Equation  24 
with  Xe,Ye,  Ze  replaced  by  the  coordinates  of  the  centroid  of  the  fragment  at 
the  first  station  and  using  the  time  interval  between  station  flashes  for  At. 

For  those  images  where  a  match  can  be  made  at  each  station  and  between 
stations,  this  method  is  more  accurate  than  the  first.  However,  it  is  more 
difficult  to  implement  and  takes  a  much  longer  time  to  evaluate  since  more 
calculations  are  involved. 


III.  AUTOMATING  THE  FILM  DIGITIZATION  PROCESS 


The  paramount  shortcoming  in  the  analysis  of  film  is  the  tedious  work 
which  still  must  be  performed  by  hand.  The  person  analyzing  the  film  becomes 
fatigued  after  a  short  period  of  time  so  that  many  of  the  images  are  ignored 
and  mistakes  are  more  readily  made.  Consequently,  an  automatic  technique  for 
performing  the  digitization  would  lead  to  a  savings  in  time  for  the  analyst 
and  a  substantial  improvement  in  accuracy. 

A  proposed  method  for  automating  the  film  digitization  process  is  based 
on  utilizing  a  video  camera  to  scan  the  film  in  a  systematic  fashion.  Because 
the  angular  camera  view  is  restricted,  a  mechanical  device  is  required  which 
can  move  the  film  so  that  adjacent  sections  of  the  film  are  positioned  into 
the  field  of  view  of  the  camera  in  a  known  sequence.  Preferably,  the  film 
sections  do  not  overlap  and  all  of  the  film  is  processed.  The  video  camera 
records  the  light  intensity  level  (gray  level)  of  the  backlighted  film  which 
is  digitized  using  an  analog  to  digital  (A/D)  converter.  The  digitized  signal 
is  then  transferred  to  a  computer  for  analysis. 

The  software  (computer  program)  for  controlling  the  analysis  should 
include  a  routine  which  will  determine  which  dots  (pixels)  are  to  be  asso¬ 
ciated  with  each  image.  For  each  image  which  is  resolved,  the  area  and  the 
five  points  describing  the  image  are  determined  analogous  to  the  method  des¬ 
cribed  in  Section  II.  Incorporated  in  the  software  package  should  be  the 
software  which  has  already  been  developed  for  determining  the  fragment  mass 
and  the  velocity. 

This  method  of  analyzing  the  X-ray  film  has  an  added  benefit.  Since  the 
equipment  that  processes  the  signal  from  the  video  camera  is  recording  inten¬ 
sity  levels  (gray  levels),  the  variation  in  the  intensity  is  a  measure  of  how 
much  metal  attenuated  the  X-ray  radiation  from  the  pulsed  tubes.  That  is,  the 
thicker  the  fragment,  the  more  X-rays  are  blocked  and  the  image  on  the  film 
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reflects  that  attenuation.  Therefore, 


the  mass  can  be  estimated  directly. 


3,4 


IV.  THE  AUTOMATIC  DATA  ACQUISITION  SYSTEM  (ADAS) 


A  block  diagram  of  the  ADAS  is  shown  in  Figure  6.  The  X-ray  side  of  the 
system  (the  X-ray  tubes)  is  the  same  as  current  conventional  range  setups.  The 
tubes  are  pulsed  using  the  same  triggering  mechanism  as  now  and  with  the  same 
equipment.  What  is  different  from  the  conventional  setup  is  that  the  film  has 
been  replaced  by  a  set  of  video  cameras.  The  intensifier  screens  (also  used 
with  the  conventional  film  system),  which  are  located  in  front  of  the  cameras, 
convert  the  X-ray  radiation  that  impinges  on  the  screen  to  visible  light.  The 
images  of  the  fragments  appear  as  dark  areas  on  the  intensifier  screen  since 
the  fragments  attenuate  or  block  the  radiation.  A  lens  focuses  the  view  of  the 
intensifier  screen  onto  the  photo-cathode  of  a  micro-channel  plate  (MCP)  image 
intensifier  tube  when  the  MCP  is  turned  on  (triggered  by  the  same  pulse  that 
triggered  the  associated  X-ray  tube)  by  the  operation  of  a  control  unit.  In 
addition  to  acting  as  a  shutter,  the  MCP  amplifies  the  light  focused  on  its 
photo-cathode.  The  gain  of  the  MCP  amplifier  is  approximately  1000  which  is 
sufficient  to  compensate  for  the  extremely  short  exposure  time  that  the  MCP  is 
gated  on.  (The  exposure  time  is  comparable  in  length  to  the  duration  of  the  X- 
ray  pulse  of  radiation  which  is  approximately  70  nanoseconds  and  precludes 
blurred  images  of  the  fragments.)  The  amplified  light  signal  is  coupled 
through  a  fiber  optic  cable  to  the  vidicon  tube  of  the  video  camera. 

Figure  7  is  a  schematic  representation  of  a  vidicon  camera  tube.  The 
signal  electrode  is  operated  at  a  positive  voltage  with  respect  to  the  photo- 
conductive  target  (PCT)  which  operates  at  the  cathode  (near  zero)  voltage. 
Pefore  the  MCP  has  been  gated  on,  the  back  side  of  the  PCT  is  charged  by  the 
action  of  a  scanning  electron  beam  which  is  supervised  by  the  control  unit. 
When  the  MCP  is  gated  on,  the  scanning  electron  beam  is  blanked  off  and  the 
light  pattern  projected  by  the  optic  fiber  cable  onto  the  PCT  causes  its 
conductivity  to  increase  in  the  illuminated  areas.  This  increased  conductivity 
causes  the  adjacent  areas  on  the  rear  of  the  PCT  to  charge  to  more  positive 
values.  The  MCP  is  then  gated  off  so  that  X-ray  radiation  from  any  tubes  that 
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Figure  7  Schematic  Representation  Of  A  Vidicon  Tube. 


are  subsequently  pulsed  will  not  affect  the  image  pattern  stored  on  the  PCT. 
To  process  the  image  pattern,  the  scanning  electron  beam  is  turned  on.  The 
signal  resulting  from  th<=>  electron  beam  interacting  with  the  positively 
charged  areas  of  the  PCT  is  capacitively  coupled  to  the  signal  electrode 
(designated  signal  output  in  Figure  7).  This  signal  after  amplification  is  the 
output  signal  from  the  video  camera. 

The  output  signal  from  the  video  camera  is  coupled  to  the  A/D  converter. 
The  A/D  converter  periodically  samples  the  output  signal  from  the  video  camera 
and  converts  the  amplitude  of  the  signal  at  that  instant  to  a  digital  value. 
This  value  corresponds  to  the  intensity  of  the  light  at  that  position  in  the 
image  pattern  on  the  PCT.  This  action  continues  until  two  complete  scans  (one 
frame  for  standard  interlaced  video)  of  the  PCT  have  taken  place.  The  digi¬ 
tized  signals  are  stored  in  contigous  memory  locations  of  the  frame  store 
memories  (see  Figure  6)  as  one  frame.  The  memory  can  then  be  dumped  under 
control  of  the  minicomputer  to  a  mass  storage  device  such  as  magnetic  disk  for 
"permanent"  storage.  The  contents  of  each  frame  stored  in  memory  can  also  be 
transferred  to  a  digital  to  analog  (D/A)  converter  which  will  reconstruct  the 
image  pattern  that  appeared  on  the  MCP.  The  output  of  the  D/A  converter  can  be 
displayed  on  a  video  monitor  permitting  visual  examination  of  the  fragment 
image  pattern. 


The  digitized  data  contained  in  the  frame  store  memories  (or  on  disk) 
recorded  by  each  video  camera  can  be  analyzed  using  the  methods  described  in 
Section  II  and  III.  Since  the  image  pattern  in  each  frame  represents  a  single 
"snapshot"  of  the  action  and  since  the  data  is  already  in  the  correct  form  to 
be  used  by  a  computer,  the  analysis  will  be  much  easier  than  that  required 
when  film  is  used  to  record  the  fragment  images. 

The  resolution  of  the  ADAS  is  limited  primarily  by  the  MCP  image 
intensifier  tube.  It  is  also  dependent  on  the  size  of  the  area  to  be  viewed. 
For  a  500  x  500  millimeter  viewing  area,  the  resolution  should  be  about  1  mm 
at  the  bright  Darts  of  the  image.  Technological  advances  will  likely  imnrove 
the  resolution.  If  the  resolution  is  not  deemed  adequate,  two  advantages  to 
the  ADAS  should  be  taken  into  consideration.  The  first  advantage  is  that  the 
image  data  is  handled  entirely  electronically,  thereby  eliminating  human 
error.  The  second  advantage  is  that  single  "snapshot"  views  of  the  action  are 
obtained  -  eliminating  errors  due  to  overlapping  image  regions  which  are  often 
a  problem  with  multiflash  sequences  on  flash  radiographs. 
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V.  SUMMARY  AND  CONCLUSIONS 


This  report  presents  three  different  levels  of  analytic  techniques  for 
determining  mass  and  velocity  of  dynamic  fragments  (fragments  in  motion).  The 
first  level  describes  recent  advancements  for  analyzing  images  on  flash  radio¬ 
graphs.  This  advancement  has  reduced  the  amount  of  physical  labor  and  inaccu¬ 
racies  but  the  procedure  can  still  be  made  less  tedious  by  introducing  the 
system  which  constitutes  the  second  level.  This  second  level  consists  of  a 
proposed  technique  for  eliminating  the  need  to  digitize  the  film  manually. 
This  advancement  is  entirely  feasible  but  reouires  the  procurement  of  some 
hardware  and  the  development  of  some  computer  software  for  the  necessary 
manipulation  of  the  digitized  data.  The  third  level  is  the  development  of  the 
ADAS  which  is  the  ultimate  in  automatic  recording  of  dynamic  fragment  images. 
Tnis  system  will  minimize  the  tedious  handwork  for  the  analyst  and  will  signif¬ 
icantly  increase  the  accuracy  of  the  results.  The  implementation  of  the  ADAS 
will  require  procurement  of  hardware  and  the  development  of  additional  compu¬ 
ter  software.  The  installation  of  the  video  cameras  will  require  modifications 
to  existing  firing  ranges  since  the  cameras  will  require  more  space  than  the 
film  cassettes.  It  may  be  found  that  the  ADAS  will  be  useful  for  studying 
other  physical  phenomena  such  as  shaped  charge  jet  and  self-forging  fragment 
forma ti on . 
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